#!/bin/bash
# Start the scamp program
./scamp init.cfg run.cfg
# Generate the R mass balance output 
cat plot_run.R | R --vanilla
# Generate the output for each time step
TIMESTEPS="00050 00100 00150 00200 00250"
#TIMESTEPS="00200 00400 00600 00800 01000"
for timestep in $TIMESTEPS ; do
    cat plot_timestep.R | sed -e s/"00000"/"${timestep}"/ | R --vanilla
done
